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Abstract 

We study the supersymmetric hybrid model near and after the end of infla- 
tion. As usual, we reduce the model to a purely scalar hybrid model on the 
level of the classical fields. But on the level of quantum fluctuations and their 
backreaction we take into account all superpartners of the waterfall field in 
a large-A approximation. The evolution after slow roll displays two phases 
with a different characteristic behaviour of the classical and fluctuation fields. 
We find that the fluctuations of the pseudoscalar superpartner are of partic- 
ular importance in the late time phase. The motion of the waterfall field 
towards its classical expectation value is found to be very slow and suggests 
a rather flat potential and a stochastic force. 
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1 Introduction 



The hybrid model has been proposed by Linde [H [21 [3] as a possible infla- 
tionary scenario, several variants of the model have been discussed recently 
[H El El [71 El [HI [in] ■ Here we consider the model in the context of preheat- 
ing after inflation, and not inflation itself. This period has many interesting 
aspects and has received a wide attention [5l[6l[7t[8l[9l[T0l[TTl[T2l[T3l[Til 
[15l|16l|l7l[l8l[l9l[20l|2ll|22]. The hybrid model may arise naturally in the 
context of supersymmetry and supergravity [23l [H Ull |13l [HI [3 [9]. Another 
aspect of the type of model investigated here may be that it simulates a 
second order phase transition, where the change in the effective potential 
does not arise from a decrease in temperature but is mediated by an effective 
field. It thereby replaces models [2ll [251 [13 [IB] where a rapid decrease of 
the temperature is simulated by an instantaneous quench. 

The supersymmetric hybrid model [23l [1] is often reduced to a purely 
scalar hybrid model. This is generally done on the level of classical fields. 
Loop effects are invoked sometimes in order to generate an effective potential 
for the inflaton field, replacing an explicit mass term. In the period right 
after inflation, at the time when the phase transition of the waterfall field 
takes place, and the subsequent preheating period, quantum fluctuations are 
produced abundantly, by spinodal decomposition or by parametric resonance, 
and their backreaction has to be taken into account. In previous publications 
p!3l [TOl [26j the quantum fluctuations of both the inflaton and the waterfall 
field and their backreaction on the classical fields were discussed in different 
approximations. Here we would like to investigate the role of the fluctuations 
of the superpartners, as appropriate for a supersymmetric model. The first 
studies of a supersymmetric quantum field theory out-of-equilibrium have 
appeared some years ago [271 [28], but they have not been put into the context 
of cosmology. Here we consider the supersymmetric hybrid model which is 
one of the standard multifield models of inflation. 

In a model with a spinodal instability it is essential to include the backre- 
action of the modes onto themselves. It is this backreaction which stabilizes 
the system by shifting the negative squared masses back to positive val- 
ues. Such a backreaction can be included in the most simple form either 
in the large-A^ or the Hartree approximations. The Hartree approximation 
would be much more involved; also, it is not clear how to treat the fermion 
fields consistently. Here we use, as a first approach to the problem, a large- 
N approximation, along the lines of [271 (Mj- The original model contains 



2 



an f/(l) — S0(2) symmetry for the waterfall field, which we elevate to an 
SO{N) symmetry as a basis to the large- iV approximation. 

The plan of the paper is as follows. In section [2] we present the basic 
model, its large-A^ version and the dynamical equations in unrenormalized 
form. Renormalization is discussed in Appendix |Al The numerical approach 
and the choice of parameters are discussed in section [31 The results of the 
numerical simulations are presented and discussed in section HI We end with 
a summary and conclusions in section [5l 

2 The model and basic dynamical equations 

The supersymmetric hybrid model [231 lH] is usually based on the superpo- 
tential 

Ty(^,$i,$2) = /«5(<l>i<l>2-/iV'^) • (1) 

By various arguments the model is then reduced to the standard hybrid 
model involving two fields: the infiaton field (p, as the scalar field of the 
supermultiplet S, and the "waterfall" field x ^.s the remnant of the two 
superfields $i and $2- The Lagrangean takes the form of the standard hybrid 
model 

C = + ^d.xd^x - ^mV' - ^gVx' - ^(x' - v'f . (2) 

The only aspect of this Lagrangean that links it to the supersymmetric model 
is the relation between the couplings A and g: 

A = ^ = (3) 

The vacuum expectation value v is given by f = /i/v^- Furthermore the 
mass m is often taken to be zero and the potential for the field is generated 
by the quantum fiuctuations of the waterfall supermultiplet. In Ref. [26] A. 
Heinen and one of us (J.B.) have considered this purely scalar hybrid model 
including the one-loop quantum backreaction in the Hartree approximation. 
However: even if we reduce the supersymmetric model to the scalar hybrid 
model on the level of classical fields, there are more quantum fiuctuations 
than just those of the two scalar fields, namely the quantum fiuctuations of 
the pseudoscalar and fermionic superpartners. These are dismissed in the 
usual treatment of the model, and it their impact that we want to study. 
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When taking into account the quantum fluctuations in a theory with 
spontaneous symmetry breaking one has to go beyond the one-loop approx- 
imation: the squared masses of the fluctuations can become negative, and 
this leads to a fast breakdown of the system due to an exponential increase 
of the fluctuations. In the large- iV and Hartree aproximations the backreac- 
tion of these modes onto themselves stabilizes the evolution, and this seems 
to be a sound feature of these approximations. The Hartree approximation, 
when applied to all flelds becomes quite involved, and there are conceptual 
problems with incorporating the fermion flelds. 

An approximation that has already been formulated for a supersymmetric 
system is the large- A^ approximation [27|. Here we will consider this 
approximation, but the original model with the superpotential ([T]) does not 
display any large A^ in which we could expand. But of course there may be 
generalizations like the one discussed in Ref. [129] where the inflaton couples 
to a fleld with a higher symmetry group. Speciflcally, we will consider here 
a large- A^ extension of the original model, based on the potential 

K / ^ N \ 

W{%, $1, . . . , = -^$0 (^E '^l - y^' j ■ (4) 
The previous model is recovered for A^ = 2 by the identiflcation 

$1 ^ -1=($^ +^$2) = $ , (5) 

<^2 ^ -^($1-2<I>2) (6) 

The latter notation is used in the work of Dvali et al. |23j. Note that all flelds 
are complex superflelds. The original model has a [/(I) ~ 5*0(2) symmetry 
$ — > exp(i7)$, $ — > exp(— Z7)$ which here is converted into an SO{N) 
symmetry. 

The auxiliary flelds are given by 

dW K 



The scalar potential becomes 



The Yukawa part of the fermion Lagrangian is given by 



N 

V'o [M'^ - 75) + + 75)] 

k=l 

f; K(i - 75) + viii + 75)] ^0 (10) 

k=l 
N 

^ ^fe bo(l - 75) + <^o(l + 75)] V'fe ■ 
vA^ fe=i 

The large- iV hmit is obtained by introdTicing two spatially homogeneous 
classical fields (po and 0i, and N + 1 fiuctuation fields rji via 
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The classical fields are taken to be real. 

The large- part of the bosonic Lagrangian is given by 



N 



k=2 
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,2/2 



K, 
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V ] — K 



N 



k=2 



(12) 



while the fermionic part becomes 

N 

C. . „ = — 



1 1^1- 

fe=2 



(13) 



Note that in leading order of large- we only consider the N —\ "transversal" 
quantum fluctuations r/f^.i/jk. k = 2..N. 
We now introduce the real fields ak,bk. 



rjk ^ -^{ttk + ibk) k^2..N. 



(14) 
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The part quadratic in the bosonic quantum fluctuations becomes 
1 ^ 

^ E [di^akd^ak + d^bkd^k - «'0o(4 + bl) 



1 ^ 

'NjAuct 



2. =2 



(15) 



So on the tree level the masses are given by 



ml = 2k'<I>1 + ^{cI>1-v'), (16) 

= 2k (18) 
satisfying the supersymmetry relation 

^{-l)^(m^^O. (19) 
We now introduce a Gaussian wave functional, so that 

< ai{x)aj{y) > = < a{x)a{y) > , (20) 

< bi{x)bj{y) > = < 6(x)6(y) > Sij , (21) 

< a,(x)6,(y) > = 0, (22) 
< iJi{x)i^j{y) > = < ^{x)7jj{y) > Sij (23) 

and all the higher correlation functions reduce to the two-point ones. We 
now can replace 

E < OfeOfe >= (AT - 1) < >~ AT < > (24) 

k 

and similarly for all terms which are second order in the fluctuations. The 
quartic term needs some care: we have 

< E(«fc + 2?afe&fe - bl){al + 2iaibi - bf) >= N^{< > - <b^ >f + 0{N) 

k,l 

(25) 

as for correlations between some Ofe and some a; we get a Kronecker 5ki and 
for such correlations only one summation remains, yielding a factor N instead 
of N^. 
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Using the Gaussian factorization and treating all transverse fluctuations 
as dynamically identical we obtain the bosonic part of the Lagrangian: 

-^(a^-ft^f, (26) 
while the fermion part takes the form 

-^^N,ferm = (^1^,^!^ " ^\f^4>^ ■ (27) 

The effective masses now become 

Ml = 2K'^(t)l + k\(1)1-v^) + [< > - Kb"^ >) , (28) 

Ml = 2K^(t)l - ^^^{<^\ - ^'^) -n^{<a^ > - <b^ >) , (29) 

Mj = 2k^(I)1, (30) 

or 

Ml = M} + M^, (31) 

Ml = Mj-Ml. (32) 

The supersymmetry sum rule for the masses is still satisfied. 
The mass M'^_ satisfies a gap equation 

Ml_ = K^{(l>i - v^) + K-'(<a^>-<b^>) , (33) 

where the right hand side is a functional of 0i and M'^_. The gap equation 
has to be solved at f = 0, later on it is then satisfied automatically. We define 

m% = M'j{0) J = a,b,f (34) 

and the potentials 

V,{t) = M,{t) - m% . (35) 
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The equations of motion for the classical fields are obtained as 

00 + 2/^^001 + 2kVo(< > + < >) + i« < '0V' > = 0, (36) 

v2 

01 + 2^20^01 + «;Vi(0? - -y^) + k20^(< a2 ^ _< ^2 = Q. (37) 

The second equation may also be written as 

0i + M^0^ = O. (38) 

Of course it remains a nonlinear equation and the equations for 0o and 0i 
form a coupled system. 
The fluctuations satisfy 

'd-Aa + Mla = 0, (39) 
b-Ab + Mlb = 0, (40) 
{i^d^ - V2K(Po)ip = 0. (41) 

Unlike the case of the hybrid model in the Hartree approximation the system 
of equations for the different fluctuations are uncoupled in the large- A?" limit. 
This would change if 0i were allowed to be complex. 
We expand the fluctuations into mode functions as 

«(x,t) = /^(a(k)/,"(t)e^'^- + at(k)/r(t)e-'^-), (42) 

bM = /^(6(k)/,^(t)e^'^- + 6t(k)/,^*(t)e-^-), (43) 

^ [c(k, s)U,{k, t) + ct(-k, s)K(k, s, t)] e^"^-^ . (44) 

The mode functions fl^^{t) satisfy 

'if + k^ff + Ml,jf = (45) 
and have the initial conditions 

ff{0) = 1, (46) 
ff{0) = -icua/w, (47) 
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with Ujo = ^k'^ + m^Q, j = a, b, f. 
We write the spinors Us(k,t) as 



f/,(k,t) =iVo idt + Hi^it) ft{t) 



K(k,t) =iVo idt + n^^{t) gt{t) 



Xs 




Xs 



with the Fourier-transformed Hamiltonian 

7ik(t) = cxk + Mf{t)(3 . 

For the two-spinors Xs we use hehcity eigenstates, i.e., 

k(TX± = ±X± ■ 

The mode functions fk{t) and gtii) only depend on k 
second order differential equations 

tMfit) + e + M){t) 



(48) 
(49) 

(50) 
(51) 

|k|; they obey the 



f^{k,t) 

gtit) 








The initial conditions are 

ftio) 
gtio) 

so that gtit) = ftit). 



(52) 
(53) 

(54) 
(55) 



Then the mode sums or fluctuation integrals take the form 

d^k 



< a" > 



<b'> 



< tp-ip > 



-2 



(27r) 32^,0 
d^k 

(27r)32cu;,o 

d^k 



\fM' 



(27r)32a;/o 



2uJfo - 



2k2 



ujfo + m/0 



\ftit)\' 



(56) 
(57) 
(58) 



All the fluctuation integrals are divergent. The procedure to separate 
them into finite integrals over the mode functions and renormalized finite 
parts of the divergent leading order contributions has been described in 
Refs.|30[ [31]. The renormalization is discussed in Appendix lAl 



9 



3 Numerical procedure and choice of param- 
eters 



We have implemented the equations of motion derived in the last section into 
a numerical code. After separation of the divergent parts and renormalization 
the computation of the fluctuation integrals reduces to finite integrals and is 
straightforward. The procedure is described in detail in Refs. [301 EI] • The 
choice of the momentum grid and of the time steps depends of course on the 
parameters of the model and on the initial conditions. For our choices (see 
below) we have typically extended the momentum integration up to Pmax = 15 
and used a grid of 700 momenta, slightly concentrated towards p = 0. The 
reliability was checked using energy conservation and by the constancy of 
the Wronskians. With small values of k the time evolution is very slow. 
Nevertheless the integration of the mode functions at high momenta requires 
sufficiently small time steps. We have chosen At ~ 10"'* and followed the 
evolution up to times of tmax — 60000. 

The original model has two parameters, k and /i. We have chosen /i = 1, 
which then defines the mass and inverse length scale. Then v = ^l.! \fK = 
l/y/n. This fixes the units we have chosen for the plots. The renormalization 
scale was taken to he m = kv = /i-^K. For each run we have to specify 
two more parameters, 0o(O) and 0i(O). We have to give the waterfall field 
a small initial value because otherwise it remains zero forever. The value 
cannot be chosen "infinitesimally small" like 10~^, as in the absence of Hubble 
expansion and of an explicit infiaton mass, the system remains at 0o = 0o(O) 
for a very long time. The fiuctuation integrals are nonzero at t = and can 
in principle initiate a time evolution of (po, but for small initial excitations 
their numerical values are too small and for larger ones there is a cancellation 
between bosonic and fermionic contributions to the equation of motion of 0o- 
Typically we have chosen |0i(O)| ~ 0.01. 

As to the parameter k a value of k = 0.01 has been suggested in Ref. 
[23], much smaller values have been discussed in Ref. [32], the simulations in 
Ref. [13] were presented for k = 0.001. In any case k is a small parameter, 
and we have performed simulations with k = 0.1, 0.01 and 0.001. As the 
time evolution gets slower with decreasing values of n we present most of 
our results for k = 0.1. The qualitative features remain the same for smaller 
values, but a detailed study would require very long CPU times. 
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4 Numerical results 



4.1 Time evolution I: the slow roll period. 

This period takes a time ti which ranges between 50 for n = 0.1 and low 
excitations and several thousand for k = 0.001 and large initial value of 0o- 
The evolution is mostly classical, the fluctuations remain very small. This 
period ends once the mass squared of the scalar (a field) fluctuations gets 
negative and the system enters the spinodal regime; this is the usual end of 
the slow roll phase, marked in the figures as a vertical line at t = 402 for 
0o(O) = 4 and at t = 1374 for 0o(O) = 64. The behaviour of the masses in 
this transition region is displayed in Figs. [1] and [2] for small and large initial 
exctitations, respectively. gets negative first, then the strong increase in 
the a field fluctuations drives Aif to negative values. Subsequently A^^ and 
A4l oscillate around zero with alternating signs. After a few oscillations Ail 
remains positive while Ail continues between positive and negative values. 
We have entered the second stage of evolution. 

The point where Ai^ becomes negative marks the end of inflation. At 
this point the energy density S of the classical fields is transferred to the 
quantum fluctuations, and the pressure, which was essentially equal to —£ 
in the inflationary phase becomes the one of a massive or massless gas (see e.g. 
[33] of an out-of-equilibrium analysis). This transfer of energy density can 
be seen in Figs. [3] and H] for small and large initial excitiations, respectively, 
while it takes a very long time for small excitations. We will see, however, 
that the evolution of the classical fields towards the classical minimum of the 
potential will take a long time and is not as instantaneous as often assumed. 

If the Hubble expansion were taken into account, the evolution during 
the slow roll stage would be modified essentially. So our numerical results 
for this period are not really relevant, they mostly set an initial condition for 
the second stage, which is characterized by the emergence of the quantum 
fiuctuations. 

4.2 Time evolution II: the intermediate period. 

In the intermediate period there are two qualitatively different kinds of evolu- 
tion. For small excitations 0o(O) the waterfall field 0i decides to move either 
into the positive or the negative direction and attains some average value 0i 
whose absolute value is smaller than the tree level vacuum expectation value 
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V. We call this the "broken symmetry phase". For large excitations 0o the 
waterfall field oscillates around zero. We call this the "symmetric phase". 
Such a phenomenon has similarly observed in Ref. [26], we discuss it in more 
detail below. We display the behaviour of the classical fields during this 
second and the early third period in Figs. [5] and [6] for low and high initial 
excitation, respectively. 

The quantum fiuctuations start developing right after entering the spin- 
odal regime. The fermion fiuctuations remain very small throughout. The a 
amplitudes have initially developed exponentially in a low momentum band 
and J-'a has reached large values right away. Subsequently these fiuctuations 
develop only very slowly, though Ai^ oscillates around zero and therefore be- 
comes negative periodically. This entails an exponential evolution in the low 
momentum band. However, the periods of exponential evolution alternate 
with periods where Ai"^ > 0; these intermittent periods modify the initial 
conditions for the next exponential evolution. So we have neither parametric 
resonance nor a spinodal evolution and the fiuctuation integral J-'a does not 
evolve significantly. 

The squared mass of the b field fiuctuations oscillates around some pos- 
itive value. So they could develop by parametric resonance. A resonance 
peak can indeed be seen at the beginning of, and during the second period. 
However, the amplitude of the mass oscillations is small and is of course not 
cleanly periodic. So the resonance does not evolve efficiently. 

We display the evolution of the fiuctuation integrals during the second 
and early third period in Figs. [7] and [8], for small and large initial values of 
00, respectively. 

Right after the intermediate period the structure of the momentum spec- 
tra still displays the same features as in the early stage of this period. They 
just have become somewhat irregular. They are displayed in Figs. [TTlandfT^ 
for small and large initial values of 0o- 

4.3 Time evolution III: The late time behaviour 

The intermediate period ends with a second strong evolution of the quantum 
fiuctuations, this time including those of the b fields. With the decrease 
of the field (^o, the increase of 0i and of < > the squared mass Aif is, 
at some point, driven again to negative values, see Eq. fl29|) . The time at 
which this occurs is marked by a vertical line, at t = 1733 for 0o(O) = 4 and 
t = 2438 for 0o(O) = 64. Of course the increase of the b fiuctuations drives 
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the b mass back to positive values, but then the a mass gets negative again, 
and so on. This behaviour at the onset of the late time regime is displayed in 
Figs. [131 The vertical line marks the time where Aif starts taking negative 
values again. 

Now the b field fiuctuations develop significantly. This entails strong fiuc- 
tuations in all squared masses and a characteristic change in the evolution 
of the waterfall field. While in the intermediate period it keeps oscillating 
more or less regularly around zero (symmetric phase, high excitations) or 
some finite value (broken symmetry phase, low excitations), it now starts 
shifting towards the classical minimum (pi = v in quite an irregular motion. 
The shift is very slow and the motion does not display any periodic oscilla- 
tions. So the out-of-equilibrium effective potential (if at all one may use such 
a term) seems to be quite fiat, and the motion seems to be determined by 
stochastic forces instead of well-defined harmonic forces. This is displayed in 
Figs. [H] and [151 The stochastic behaviour may originate on the one hand 
from a parametric resonance that is strongly disturbed by the presence and 
variation of different time scales. On the other hand the oscillations of A^^ 
and Ail around zero naturally lead to a diffusion process by the alternation 
of oscillating and exponential time evolution. 

In this late time period both Ai^ and Ail oscillate around zero. So, 
for the low momentum modes the effective frequency takes real and imagi- 
nary values in intermittent time intervals. As before for the fluctuations of 
the a field we get now for both scalar fields a process where an exponen- 
tial behaviour alternates with an oscillating evolution, which changes initial 
conditions for the next exponential evolution. This time the amplitudes of 
oscillation for Atl/f, are larger, the process becomes more effective and leads 
to a strong increase of the fluctuation integrals, as mentioned above. At the 
same time the momentum spectra loose entirely their characteristic features 
related to parametric resonance and/or spinodal evolution. They are concen- 
trated at small momenta and fall off rapidly with momentum. These spectra 
are displayed in Fig. [161 for 0o(O) = 4 and in Fig. [TTlfor 0o(O) = 64. 

This third period of evolution shows novel features that deserve to be 
analyzed in more detail. Without the pseudoscalar superpartner, the system 
remains in the "symmetric" and "broken symmetry" phases |26j, here the 
waterfall field always evolves towards of one of its classical minima at late 
times. The stochastic behaviour contrasts with the spinodal or parametric 
resonance regimes and seems to represent a new kind of behaviour specific 
of multifield quantum systems. 



13 



4.4 A phase transition 

As we have discussed previously the system behaves in two different ways 
in the intermediate stage: for small excitations the field 0i shifts towards 
the direction of one of the classical minima and oscillates in an irregular 
way around some mean value which depends on the initial excitation. With 
increasing initial excitation this mean value becomes smaller and above some 
critical value it becomes zero. This can be seen as a phase transition, with a 
symmetric phase for large excitation, i.e. large energy densities. We plot, in 

Fig. m 

this average value as a function of 0o(O), for k, = 0.1. As the intermediate 
period only lasts for a finite time, this average is not really well-defined and is 
read off by inspection. This explains why the values are somewhat scattered. 
The plot looks like the phase diagram of a second order phase transition, with 
00 (0) as the "temperature" and (pi as the order parameter. The diagram may 
be compared with similar figures in other work on models with spontaneous 
symmetry breaking, like Fig. 11 of Ref. [26j, Figs. 6-8 in Ref. [34j and Fig. 
4 in |35J. In contrast to the models studied previously, for the present model 
the waterfall field always ends up at |0i| ~ f as t ^ oo, but this is a process 
that takes a long time. 

4.5 The effective potential 

As we have noticed, the motion of the waterfall fields towards its tree level 
expectation value does not seem to be driven by any strong force. The slow 
and somewhat irregular motion rather suggests an effective potential that is 
quite flat, with a small inclination towards 0i = ±f . Indeed it is well-known 
that for 0^ theory with SO{N) symmetry and a Mexican hat potential the 
central region |0| < f becomes fiat, once quantum corrections are taken into 
account in the large- limit. Here a similar phenomenon takes place. The 
effective potential for our model in the large- A^ limit is derived in Appendix 
|B]and is displayed in Fig. [13 

The region in the center is fiat in the 0i-direction, it is the region B 
defined in Appendix [B] where the gap equation would yield the unphysical 
solution Ail < 0. At 00 = it is bounded by 0i = ±v, i.e. by ±1/-\/k in 
our mass units. One sees that the spinodal region has disappeared and the 
effective potential has become fiat in the central region. Of course this is the 
equilibrium effective potential, here its flatness seems to be "realized dynam- 
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ically" by the motion of the waterfall field, as it was found in a somewhat 
different way for the SO{N) model in Ref. pS] . 

5 Summary and Conclusions 

We have studied here a supersymmetric hybrid model in a large- iV approx- 
imation. The motivation for this research was the interest in the role of 
quantum fluctuations in the evolution near and after the end of the slow 
roll regime. In generalizing previous work we have taken into account the 
fluctuations of all superpartners of the waterfall field. While the fermion 
field fluctuations remain small in general and they mainly play their role 
in the renormalization, the fluctuations of the scalar and pseudoscalar fields 
show an interesting interplay, resulting in a peculiar behaviour of the clas- 
sical waterfall field which is different from the behaviour found in previous 
analyses: 

(i) the evolution after the end of slow roll displays two stages: in the first 
(intermediate) stage the behaviour is similar to the one found previously for 
the hybrid model with quantum backreaction [26], only the scalar field fluctu- 
ations develop and the system goes to a symmetric phase at high excitations, 
and to a broken symmetry phase at low excitation. In the second stage (late 
time) the fluctuations of the pseudoscalar superpartner start developing and 
the waterfall fields moves towards its classical vacuum expectation value. 

(ii) the dynamics of the fluctuations is, apart from an initial spinodal 
evolution of the scalar field fluctuations, neither dominated by the spinodal 
instability nor by parametric resonance. On the one hand a clean parametric 
resonance is suppressed by the presence of multiple time scales, there is 
no simple periodic motion in the background fields. On the other hand 
the oscillations of the effective squared masses Ai^ and Ail around zero 
lead to a kind of diffusive behaviour for the (dominant) low energy modes. 
Within short time intervals exponential and oscillatory evolution alternate, 
the fluctuations and their integrals vary in a stochastic way. This behaviour 
is transmitted to the waterfall field whose motion is an irregular drift towards 
its classical vacuum expectation value. 

We think that these findings warrant a further consideration of the hybrid 
model in its supersymmetric extension. The following issues would be of 
interest and could be the subject for further studies: 

(i) how does the behaviour change in the presence of Hubble expansion, 
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i.e. in a FRW universe? This will lead to a damping of the classical field and 
of the fiuctuations. While the Hubble expansion is often neglected in studies 
of preheating, here, right at the end of infiation, it can be expected to lead to 
significant modifications. This is also the reason for which we do not want to 
draw premature conclusions from our numerical results. The techniques for 
handling the dynamics in an expanding universe, including renormalization, 
have been developed previously 157]. 

(ii) does the peculiar dynamical behaviour lead to imprints on the CMBR 
spectrum? In particular, nongaussianity may arise in such multifield models 
from the evolution near the end of infiation, see [22] and references therein. 

(iii) is the peculiar behaviour a genuine property of a supersymmetric 
model with its specific structure of the mass terms, or is it simply a conse- 
quence of the presence of many different time scales, or both? 

We plan to study these questions in the near future. 
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A Renormalization 

The fluctuation integrals introduced in section [2] are divergent. In dimen- 
sional regularization we find 



Here 



■fin 



^ r 



-rfin 

■J' b 



with 



and [3T 



< > 
<b^> 

<ipi) > 



(Al) 



1 



167r2 
1 

167r2 



m 



2 

aO 



I --Am' 



Ml , m 
In 



2Mf + 4:M 



167r2 

167r2 
1 



■sub 



In 



■sub 
b 



60 



87r2 



In 



(2M) + 72^00 



I -rsub 



■rsub 

a/b 



it) 



{2TT)^2uJa/bO 



\ffm' + 



(A2) 
(A3) 
(A4) 

(A5) 



(2^ 



m 



LOfo 
2/ 



Mfit) , ■MK^)(-^/(^)-^/o 



4^/0 



2uj% 



(A6) 



We introduce the renormalization constants Zj as follows: 



Po 
Pi 

K 

V 



^000 

Zi(f)i 
Ziv 



(A7) 



The fluctuation fields a, b, ip are likewise multiplied by Zi as they belong to 
the same superfield as 0i, more precisely they belong to the same SO{N) mul- 
tiplet of superfields as 0i. The vacuum expectation value has been rescaled 
as 01 in order to preserve its tree level interpretation. 
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We first discuss the gap equation. It takes the form 



with T'^ = JF^" — J^h^- This can be written as 



fin 
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(A9) 



In order to obtain the finite equations in the same form as the unrenormalized 
ones we have to put 1 + 2K^Z'^Zll_^{rn?) equal to Z'^Zf or 



zlzl 



1 - 2fi:2J_3(m2) 
and obtain the finite gap equation 

=^2 (^2_^2_^^fin^) 



(AlO) 
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Next we consider the equations of motion. We start with the one for 0i, 
which after divison by Zi takes the form 



+ 2ZIzIk'(PI(P^ + Mt(t)i = . 



As TVI^ is already finite we have to set 



ZqZi^ 



1 . 



(A12) 



(A13) 



This entails that also the renormalized equation of motion for the field ip 
retains its bare form. We now consider the equation of motion for 0o- It 
becomes, after dividing by Zq which is equivalent of multiplying with 



^0 + Zlzl 
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We use the relations (lAip to rewrite this as 
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where — J^'^ + J^^. In order to obtain again the tree level relation with 
finite quantities we have to put 

1 + 2ZlK'^Zfl_s(m^) = ZlZ^ . (A16) 

which is the same relation we have obtained previously. 

Up to now only the products ZqZ^ and ZiZ,^ are fixed. Wc now consider 
the renormalization of the energy. The total energy density written in terms 
of the bare fields and couplings is given by 

^ = 1^1 + 1^1 + ^'<l>l<l>i + ^{<l>l-vf 

+ <la' + \k'a' + \My > (A17) 



(< a' > - < 6' >)' + (ak + (3Mf) i/j . 



4 V ^ 2 

In terms of the renormalized fields this becomes 
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+Zf < + \ea' + \My > (A18) 
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Separating the divergent parts wc have 

= -\l-M^)MU^ll. (A19) 
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Here 
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where the superscript 'sub' refers to the subtracted fluctuation integrals 
which can be found in detailed form in Refs. [311 [30] . 

The divergent parts of the three fluctuation energies combine into 
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In the total energy the fluctuations contribute a term 



(A26) 



If we want the energy density to keep its bare form in terms of finite quantities 
we have to set Zi = 1. Then in Eq. (lAlSp the second and third term have 
already their bare form in terms of renormalized quantities. The first term 
and 0Q-term in the divergent fluctuation energies combine into 



This takes the canonical form if 
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This is consistent with Eq (IA13P if 
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and this is also consistent with Eq. flAlOp if Zi = 1 as obtained previously. 
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The fourth term in the energy density combines with the seventh term 
into 



= ^Ml [0? -k" - (< o - < 6 >")' 
^'.Ml{4.l-v')-j^,Mt (A30) 

The last term in the last equation cancels with the A^i-term in the fluctu- 
ation energies. So with the choices of the Zi given above the energy density 
and the equations of motion are finite. The energy density now becomes 

^ = l4+ \4>i + + + \mi - v') - ^^Mt . (A31) 

Here the gap equation has been implemented. Indeed, variation with respect 
to yields the gap equation 

M^t) = -v' + . (A32) 

The equations of motion for the fluctuations as well as the one for 0i have ex- 
actly the same form as the bare equations if we use the masses A4l, A4l, A4^f 
and M^_. The equation of motion for 0o has the form 

^0 + 2«'0o ((pl + ) + ^J^f = . (A33) 

This equation of motion cannot be solved trivially with respect to 0o as 
explicitly contains (pQ. So the term 0o appears with a factor 



Co=l + 2— ^In^. (A34) 
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Likewise the flnite gap equation contains on both sides. At t = this is 
not a problem as it is solved by iteration anyway. But if we want to express 
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at finite t by the subtracted fluctuation integrals we have to write it in 
a modified form. We have 
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B The large-N effective potential 

The energy density for constant fields is the effective potential. This effective 
potential is obtained in the 2PI formalism by maximizing with respect to 
Ail the action as a functional of the fields and masses. In our model and in 
the large- A?" approximation this action takes the form 
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Variation with respect to the mass Ml yields the gap equation: 
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or 
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(B2) 



The solution of this equation is to be inserted into V(7V1^, (po,(pi) in order to 
obtain Kfr(0o,0i)- 

The gap equation has a solution only in a restricted part of the (po, 0i 
plane. Both Ai^ and Ail have to be positive, and such a solution does not 
exist everywhere. The boundaries of this region are traced by the conditions 
Ml = and Ml = 0. Outside this region the maximum of the variational 
potential is to be taken at the boundaries. In spontaneously broken 0^ theory 
at large N this construction leads to a flat potential in the central region 
101 < f of the "Mexican hat" (see related discussions in Ref. [38]). 

Here the boundaries of the regions where the gap equation yields the un- 
physical solutions Ml < (region A) and Ml < (region B) are determined 
by the conditions 
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In the regions A and B the masses Ml and Ml, respectively, retain 
their boundary values, i.e. zero. Then the effective potential can be given 
analytically as 
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The latter potential is independent of 0i and replaces the potential in and 
around the spinodal region. 

The conditions (dA) and (dB) do not determine the regions of validity in 
terms of the fields (po and 0i. A way of obtaining these regions in a numerical 
code is to solve the gap equation admitting negative values oi Ml and Ml, 
taking the absolute values in the logs. If Ml is found to be negative one is 
in region B, if Mf is found to be negative one is in region A. As long as 
the potential V has just one maximum as a function of M^ this recipe is 
rigorous. 

If one omits the fluctuation terms, the region where the potential V^f is 

independent of 0i is the region inside |0i| < \Jv'^ — 20o which includes the 
spinodal region. 
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Figure 1: Evolution at the transition between the slow roll and the interme 
diate period, for k = 0.1 and 0o(O) = 4. 




Figure 2: Evolution at the transition between the slow roll and the interme 
diate period, for n = 0.1 and 0o(O) = 64 . 
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Figure 3: Evolution of the classical and fluctuation energies after the slow 
roll period, for k = 0.1 and 0o(O) = 4. 




Figure 4: Evolution at the classical and fluctuation energy densities after 
slow period, for k = 0.1 and 0o(O) = 64. 
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Figure 5: Evolution of the classical fields during the intermediate period for 
K = 0.1 and 0o(O) = 4. 
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Figure 6: Evolution of the classical fields during the intermediate period, for 
K = 0.1 and 0o(O) = 64. 
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Figure 7: Evolution of the classical fields during the intermediate period for 
K = 0.1 and 0o(O) = 4. 
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Figure 8: Evolution of the classical fields during the intermediate period, for 
K = 0.1 and 0o(O) = 64. 
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Figure 9: Momentum spectra for the fluctuations of the fields a (left) and b 
(right) right after slow roll, for k = 0.1 and 0o(O) = 4. 
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Figure 10: Momentum spectra for the fluctuations of the fields a (left) and 
h (right) right after slow roll, for n = 0.1 and 0o(O) = 64. 
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Figure 11: Momentum spectra for the fluctuations of the fields a (left) and 
h (right) after the intermediate period, for k = 0.1 and 4>q{0) = 4. 
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Figure 12: Momentum spectra for the fluctuations of the flelds a (left) and 
b (right) right after the intermediate period, for k = 0.1 and 0o(O) = 64. 



33 



0.04 l—^ 1 : r 




Figure 13: Evolution of the masses at the transition between intermediate 
and late time regime, for k = 0.1 and 0o(O) = 4. 




Figure 14: Evolution of the classical fields during the intermediate period for 
K = 0.1 and 0o(O) = 4. 
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Figure 15: Evolution of the classical fields during the intermediate period, 
for K = 0.1 and 0o(O) = 64. 




Figure 16: Momentum spectra for the fiuctuations of the fields a (left) and 
b (right) at late times, for k = 0.1 and 0o(O) = 4. 
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Figure 17: Momentum spectra for the fluctuations of the fields a (left) and 
b (right) at late times, for k = 0.1 and 0o(O) = 64. 




Figure 18: Averages of the waterfall field in the intermediate period as func- 
tion of the initial value (/)o(0), for k = 0.1. 
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Figure 19: The equilibrium effective potential for k = 0.1. 
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